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Transition temperatures T c calculated using the BCS model electron-phonon in- 
teraction without any adjustable parameters agree with empirical values for quasi-2D 
cuprate superconductors. They follow from a two-dimensional gas of temperature- 
dependent Cooper pairs in chemical and thermal equilibrium with unpaired fermions 
in a boson-fermion (BF) statistical model as the Bose-Einstein condensation (BEC) 
singularity temperature is approached from above. The linear (as opposed to 
quadratic) boson dispersion relation due to the Fermi sea yields substantially higher 
T c 's with the BF model than with BCS or pure-boson BEC theories. 
PACS # 03.75Fi; 05.30.-d; 05.30.Fk; 05.30.Jp 

We provide support to a widespread conjecture (or "paradigm" ) that superconductivity in general 
is a Bose-Einstein condensation (BEC) of the charged Cooper pairs |£j (CPs) observed in magnetic 
flux quantization experiments in classical [^],[| as well as cuprate 0] superconductors. The same 
general conjecture is also often made regarding the superfluidity of liquid helium-3 in terms of 
CPs consisting of neutral-atom 3 He fermions. BEC as a statistical (as opposed to a dynamical) 
mechanism of superconductivity has been entertained, among others, by Anderson g, by T.D. Lee 
§] and by Mott [§ and their co-workers. But BEC normally occurs only for dimensions d > 2 while 
some modern superconductors are quasi-2D or even quasi-lD materials. We show, however, that 
CPs can undergo BEC for all d > 1. We further obtain reasonable critical temperatures T c without 
any adjustable parameters, thus bolstering the above mentioned conjecture even before building in 
full many-body self-consistency. As in BCS theory, fluctuations have also been neglected. More 
detailed, sophisticated treatments actually link [p|-p^[ . BEC (characterized by a bosonic condensate 
fraction) with the BCS theory (characterized by a fermionic gap), but report no attempts to 
calculate specific T c 's without adjustable parameters to compare with experiment. 

A BEC picture of superconductivity is consistent with the recent discovery of the "pseudogap" 
in the electronic density of states [fU3|-|l9|] above T c in certain cuprates, at least with one of its 
major interpretations as "pre-formed CPs" without long-range coherence or condensation, while 
in BCS theory CP formation and condensation occur simultaneously below the same T c . We here 
submit that a natural candidate for such pre-formed CPs are the nonzero-center-of-mass CPs usually 
neglected in BCS theory. 

To fix the dynamics take a 2D system of N fermions of mass m confined in a square of area 
L 2 interacting pairwise via the BCS model electron-phonon interaction Vk,k' = —V, with V > 0, 
whenever /j,(T) — hujD < (= K kf/2m), e^ 2 < fj,(T) + fiu)D, and zero otherwise, where k = |(ki — 
k 2 ) is the relative wavevector of the two particles; n(T) the ideal Fermi gas (IFG) chemical potential, 
which at T — becomes the Fermi energy Ep = Ti 2 k 2 F /2m with Uf the Fermi wavenumber; and 
ljd the Debye frequency. Striking direct evidence for significant electron-phonon coupling in high- 
temperature cuprate superconductors from angle-resolved photoemission spectroscopy (ARPES) 
experiments has recently been reported ptj| |. 

If KK. = h(ki + k 2 ) is the center-of-mass momentum (CMM) of a CP, let Ek be its total energy 
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(besides the CP rest-mass energy). The original CP (lj eigenvalue equation is then 



1 _ v y-' g(fa ~ M 9{k 2 - k F ) 
^ 2e k + h 2 K 2 /4m- E K ' 



where 6{x) is the Heaviside unit step function, and the prime on the summation sign denotes the 
conditions hip = |K/2 ± k| < (k% + k 2 ^) 1 / 2 ensuring that the pair of particles above the Fermi 
"surface" cease interacting beyond the annulus of energy thickness 2h\op> = h 2 k 2 D /m, thereby 
restricting the summation over k for a given fixed K. CPs obey BE statistics since there is an 
indefinitely large number of k values in (Q) for a given value of K. Setting Ex = 2Ep — Ax, 
a pair is bound if Ax > 0, so that (P becomes an equation for the (positive) pair binding energy 
Ax- Our Ax and A follow Cooper's notation and should not be confused with the BCS energy 
gap A(T) at T = 0. Let 17(e) be the electronic density-of-states (for each spin) in the normal (i.e., 
interactionless) iV-fermion state; in 2D it is constant, g(e) = L 2 m/2-K% = g. The Cooper equation 
(|l|) for the unknown quantity Ax can be analyzed beyond the usual zero-CMM, K = 0, case. For 
K = it becomes a single elementary integral, with the familiar JlJ solution Ao = 2hur}/(e 2 / x — 1) 
valid for all coupling A = gV > 0. For small A one gets 

Ax — > A - -hv F K + 0(K 2 ) (2) 



where vp = yf2Ep jm is the Fermi velocity. This linear dispersion relation is the 2D analog 
of the 3D result stated by Schrieffer as far back as 1964 in Ref. p2] |, p. 33 (see also Ref. p^ |, 
p. 336). Although some treatments (e.g., Ref. (^4|) of CPs more sophisticated than the original 
Cooper picture ([!]) numerically yield resonant pairs with a leading quadratic dispersion, linearly- 
dispersive resonances appear analytically from a Bethe-Salpeter equation many-body approach p5| 
to CPs in 3D — provided it is based on the BCS (where holes are treated on an equal footing 
with particles), not the IFG, ground state. In 2D, see also Refs. |^6] , p7| . It is commonly confused 
with the also linearly-dispersive sound phonons of the collective excitation sometimes referred to 
as the Anderson-Bogoliubov-Higgs mode (which for zero coupling reduces |28| to the IFG result 
nvpK I \fd). The IFG sound speed c = vp/^/d follows trivially from the zero- temperature IFG 
pressure P = n 2 [d(E/N)/dn] = 2nEpj (d + 2) via the familiar thermodynamic relation dP/dn = 
mc 2 , where E is the ground-state energy and n = N/L d — kp/2 d ~ 2 ir d / 2 d T(d/2) is the fermion- 
number density. But the simple result (||) in fact refers to actual "moving" (or "excited") CPs in 
the Fermi sea, which clearly "break up" for K > Kq as defined by Ax = 0. Both kinds of distinct 
soundwave-like solutions — moving CPs and ABH phonons — appear side by side in the many-body, 
ladder-summation scheme of Ref. |2j| . 

For Nb bosons of mass tub and energy ex = C s K s with s > and C s a constant, a BEC 
temperature singularity occurs at T c ^ for any dimension p9|j30[ | d > s in the number equation 
Nb = ^ K [e^ A_AlB )/ fcBT — 1] _1 at vanishing bosonic chemical potential fiB < when the number 
of K = bosons just ceases to be negligible upon cooling. It is given |3l| by 



T — fli 



sr(d/2)(27r) d n B 1 s/ ° 



2^/ 2 T{d/s)g d/s {\) 



(3) 



with ub = NB/L d , and g a (z) the usual Bose integrals expandable as infinite series which are 
C(<t), the Riemann zeta function of order <r, for a > 1 but diverge for a < 1. Thus T c = for 
all d < s. For s — 2 and d — 3 one has C(3/2) ~ 2.612, and since C2 = h 2 /2mB ^ then 
reduces to the familiar formula T c ~ 3.3lh 2 n 2 J 3 /rusks of "ordinary" BEC. But for bosons with 
(positive) excitation energy ex = Ao — Ax given approximately by the linear term in (|) for all K , 
meaning that s — 1 and C\ = a(d)hvp with a(d) — 2/ir and 1/2 for d = 2 and 3, respectively, the 
critical temperature T c is nonzero for all d > 1 — precisely the dimensionality range of all known 
superconductors down to the quasi-lD organo-metallic (Bechgaard) salts p^-Q. 
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The number of bosons in the boson- fermion mixture to be analyzed turns out to be both coupling- 
and temperature-dependent and it is in conserving the fermion number that a BEC-like singular- 
ity arises. As is the case for the pure boson gas, a linear rather than a quadratic dispersion 
relation is needed to obtain BEC in 2D. This emerges in a statistical model for the ideal binary 
mixture of bosons (the CPs) and unpaired (both pairable and unpairable) fermions in chemical 
equilibrium p3,p6| for which thermal pair-breaking into unpaired pairable fermions is explicitly 



allowed |37|. Assuming the BCS model interaction the total number of fermions in 2D at any T is 
N = L 2 kp/2ir = Ni + N 2 , where N% is the number of unpairable (i.e., non- interacting) fermions 
while N 2 is the number of pairable (i.e., active) ones. The unpairable fermions obey the usual 
Fermi-Dirac distribution with the IFG chemical potential /1 but the N 2 pairable ones are simply 
those in the interaction shell of energy width 2Tluid so that, if [3 = (fc^T) -1 , 

N 2 = 2 / de J\ ■ - = 2ghLo D , (4) 



which is independent of T. At hxed interfermionic coupling and temperature these N 2 fermions 
form an ideal mixture of pairable but unpaired fermions plus CPs created near the single-fermion 
energy /i(T), with binding energy Ak(T) > and total energy 

Ek(T) = 2/i(T) — Ak(T). (5) 

This generalizes the T = equation Ek = 2Ep — introduced before. 

The Helmholtz free energy F = E — TS, where E is the internal energy and S the entropy, 
of this binary "composite boson/pairable-but-unpaired-fermion system" at temperatures T < T c is 



then readily constructed 37 in terms of: a) the average number of unpaired but pairable fermions 
with fixed energy; b) Nb,k(T), the number of CPs with nonzero-CMM, < K < Kq, with the 
CP-breakup value K defined |HJ by A Ko = 0; and c) N B ^o(T), the number of CPs with zero CMM 
at temperature T. The free energy F 2 of just the N 2 pairable fermions is to be minimized subject to 
the constraint that N 2 is conserved, i.e., one seeks the minimum of F 2 — fi 2 N 2 with respect to (a), 
(b) and (c) just mentioned. The total number of pairable but unpaired fermions N 2 q(T) is then 

N 20 (T) = 2g de^ , = ^ln a , — r — r . (6) 

The relevant number equation for the pairable fermions is thus 

N 2 = N 20 (T) + 2[N Bfi (T) + n b,k{T)] = N W (T) + 2N B (T), (7) 

K>0 

where J2k">o n b,k{ t ) = X)^° >0 [e^ {B,<r(T)_2w} - is the total number of "excited" CPs (namely 
with CMM values < K < K Q ). One can rewrite E K (T) - 2pL 2 here as e K (T) - (i B (T), with 
Ek(T) = A (T) — Ak(T) > a (nonnegative) excitation energy as suggested by (||). Hole- hole 
and particle-particle CPs can be shown to have the same excitation energy ek (T) . The remaining 
unknown fx B (T) is then 

MB (r) = 2[ M2 (r)- Ai (r)] + A (T) = (8) 

for < T < T c since N Bfl {T) is negligible for all T > T c . This is precisely the BEC condition for 
a pure boson gas, although one now has a binary boson-fermion mixture. 
To determine N B (T) from (0) and (R) we use M) and see that 



^2o(T) = ^ln 

P 



I , e -/3{A (T)/2-hw D } 



I + e -/3{A (T)/2+hu D } 



(9) 
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for < T < T c . Thus 2N B {T)/N 2 = 1 - N 2Q {T)/N 2 is obtainable for < T < T c from (|) if 
A (T) were known. For this, 8(ki — kp) = 9{ek 1 — Ep) in (|l|) becomes 1 — n(£fcj where ™(£fci) = 
(e K "i + I)" 1 with £ kl = e kl - fi(T), the IFG chemical potential ^(T) in 2D being given exactly 
by fi(T) = /3- 1 lnfe^" - 1) — ► E F as T -> 0. Similar arguments hold for 6>(fc 2 - fcj?). Since for 
K = 0, fci = which implies that £k ± — £,k 2 ; (@) then provides a simple generalization to finite-T 
of the K = CP equation, namely 

1 = A I*"" d£(e- K + l)- 2 [2^ + Ao(T)]- 1 . (10) 
Jo 

Its numerical solution shows Ao(T') to decrease monotonically with T for hxed A and tkoo, and zero 
only at infinite T . (This infinite "de-pairing" temperature is obviously spurious as the BCS model 
interaction loses meaning when /i(T) turns negative at large T.) Thus also 2Nb{T)/N 2 decreases 
with T; it is plotted in Fig. 1 as 2n B (T)/n 2 , since n B (T) = N B {T)/L 2 and n 2 = N 2 /L 2 . 

Using (|J) for N 2 the fractional number of pairable fermions that are actually paired at T = 
becomes simply 

2N B (0)/N 2 = A /2hu> D = (e 2/A - l)" 1 — ► e~ 2/x (11) 

A — ^0 

for A < 2/ In 2 ~ 2.89, and = 1 (all pairable fermions paired into bosons) for A > 2.89. This fraction 
is plotted against coupling A in Fig. 1, and contrasts sharply with the "heuristic model" of Ref. pl| , 
Eq. (16), where 2Nb(0)/N 2 = 1 for all coupling. It is now more in line with BCS theory — which 
is not |}8) a BEC theory — where, in any d, a coupling-dependent fraction is estimated (Ref. |}9) 
p. 128) to be (A/huio) 2 = (sinhl/A) -2 — > 4e~ 2 / A as A — > 0. Here A (again, not to be confused 
with the CP binding energy Ao) is the T — BCS energy gap for the same BCS model interaction 
used in this Letter. It is graphed as the thin curve in Fig. 1 and is seen to be much larger than 
@ for fixed A. 
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FIG. 1. Fractional number of pairable fermions 
that are actually paired vs. coupling A for the 
present statistical model at three different temper- 
atures (thick curves) and estimated for BCS theory 
at T = as explained below (0) (thin curve). The 
number of pairable fermions with the BCS model 
interaction used is just ([lj); all of them are paired 
at T — (unrealistically) in the heuristic BEC 
model, Ref. Q Eq. (23). 
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FIG. 2. Critical BEC temperature T c in units 
of Tf for the BCS model interaction with A = 1/2 
for varying v = Uuod / n{T c ) ~ Qd/Tf for: the 
pure unbreakable-boson gas with some and with all 
fermions paired, the former being the solution of 
@ and the latter taken from Ref. 0, Eq. (17); 
for the breakable-boson gas, from Ref. J3l[], Eq. 
(18); and for the boson- fermion mixture from ( |l3| ) 
(thick full curve labeled "binary gas"). Dashed 
curve is the BCS theory T c , and cuprate experi- 
mental data are taken from Ref. pp| . 
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If the background unpaired fermions are neglected one has a pure boson gas of CPs but with 
T-dependent number density ub{T). Converting the explicit T c -formula (0) for s = 1 and d = 2 
into an implicit one by allowing ns to be T-dependent leaves 

r = ^^V^py, (12) 

since g 2 (l) = f(2) = vr 2 /6. This requires u s (T) ee N b (T)/L 2 which from (|) requires @, 
along with Aq(T) from (|l0|). Solving these i/iree coupled equations simultaneously for A = 1/2 
gives the remarkably constant value T c /Tp ~ 0.004 over the entire range of v = Twjd/Ef val- 
ues 0.03 — 0.07 typical of cuprate superconductors. On the other hand, the BCS theory formula 
t bcs _ i.ise^e-V* w ith A = 1/2 yields T C /T F = 0.005 - 0.011 over the same range of v values. 
Unfortunately, both sets of predictions are well below empirical cuprate values of T c /Tp varying 
fliofl from 0.03 — 0.09. Pure gas model results for either breakable or unbreakable bosons without 
unpaired fermions are seen in the figure to overestimate empirical T c /Tp values by factors ranging 
from two to more than two orders of magnitude. All these results are wide off the mark. 

To obtain the critical temperature without neglecting the background unpaired fermions, one 
needs the exact CP excitation energy dispersion relation Sk(T) = Aq(T) — A#(T) which is neither 
precisely linear in K nor independent of T. To determine Ak(T) we need a working equation that 
generalizes Ref. ||l[] for T > via the new CP eigenvalue equation ([n]) . For the critical temperature 
from the finite-temperature dispersion relation, besides solving for A^- (T) , one requires (^) , (Q) and 
(g). At T = T c both Nb q(T c ) ~ and /i B (T c ) ~ so that (0) leads Q to the implicit T c -equation 
for the binary gas 



V 



1 + e -{Ao(fc)/2-f}/r c 

1 + e -{A (f c )/2+ 1 /}/f c 



(1 + z/) r o(Tc) koIkj 







e [A (T c )-A K (T c )]/T c _ J 



(13) 



where quantities with tildes are in units of fi(T c ) ~ or Tp, while k = K/2\Jkp + and 
^ ee Qd/Tf- Four coupled equations must now be solved self-consistently for the exact T c for each 
A and v, including @ for A (T), and Eq. (35) of Ref. @ for both A K (f) and K (f c ). Results 
with A = 1/2 labeled "binary gas" in Fig. 2 show a huge enhancement of T c , with respect to the 
self-consistent result from (|l2|), arising from the equilibrating presence of the unpaired fermions 
and in spite of the very small number of bosons suggested by Fig. 1 for A = 1/2. 

For cuprates d ~ 2.03 has been suggested [[Il| as more realistic since it reflects inter-CuO-layer 
couplings, but our results in that case would be very close to those for d = 2 since, e.g., from (||) T c 
for s — 1 (but not for s = 2) varies little p2| with d around d = 2. Indeed, if m,Bj_ and tub are the 
boson masses perpendicular and parallel, respectively, to the cuprate planes, an "anisotropy ratio" 
itib/ 'tub _l varied from to 1 allows "tuning" d continuously from 2 to 3. 

Other boson-fermion models [||J|Jl^,^J,^,Q have been introduced, some even addressing 
d-wave interaction effects as opposed to the pure s-wave considered here, and some also focusing 
]l2] , p4] | on the pseudogap. But calculating cuprate T c values in quasi-2D without adjustable pa- 
rameters is not reported — and indeed predict T c ee in exactly 2D. 

To conclude, a statistical model treating ordinary CPs as non-interacting bosons in thermal and 
chemical equilibrium with unpaired fermions yields a boson number that rises very slowly from zero 
with coupling, and that decreases with temperature. When the CP dispersion relation is approx- 
imately linear, it exhibits a BEC of zero-CMM pairs at precisely 2D. Transition temperatures for 
the boson-fermion mixture based on the exact CP dispersion relation for the BCS model electron- 
phonon interaction are greatly enhanced over both BCS theory as well as pure-Bose-gas BEC T c 's, 
and are in rough agreement with empirical cuprate superconductor T c 's without any adjustable 
parameters. 
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